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ABSTRACT 



Context. G24.78+0.08 Al is a 20M o star surrounded by a hypercompact (HC) Hn region, driving a CO bipolar outflow, and located 
at the center of a massive rotating toroid undergoing infall towards the HC Hn region. Recent water maser observations suggest that 
the HC Hn region is expanding and accretion onto the star is halted. 

Aims. This study aims to confirm the expansion scenario proposed for the HC Hn region on the basis of recent H2O maser observations. 
Methods. We carried out continuum VLA observations at 1.3 cm and 7 mm with the A array plus Pie Town configuration to map the 
HC Hn region towards G24.78+0.08 Al. 

Results. The emission of the HC Hn region has been resolved and shows a ring shape structure. The profiles of the emission obtained 
by taking slices at different angles passing through the barycenter of the HC Hn region confirm the shell structure of the emission. The 
ratio between the inner and the outer radius of the shell, RJR a , derived fitting the normalized brightness temperature profile passing 
through the peak of the 7 mm emission, is 0.9, which indicates that the shell is thin. The deconvolved outer radius estimated from 
the fit is 590 AU. These results imply that the HC Hn region in G24 Al cannot be described in terms of a classical, homogeneous 
Hn region but is instead an ionized shell. This gives support to the model of an expanding wind-driven, ionized shell suggested by 
the kinematics and distribution of the H 2 masers associated with the HC Hn region. According to this model, the HC Hn region is 
expanding on very short times scales, 21-66 yr. 

Key words. ISM: individual objects: G24.78+0.08 Al - ISM: Hn regions - stars: circumstellar matter - stars: formation 



1. Introduction 

In the context of massive star formation, significant progress 
has been recently made with the detection of Keplerian circum- 
stellar disks around newly born B-type stars, and of massive 
toroids around young O-type stars. These results lend support 
to the accretion models that describe the formation of high-mass 
stars as a scaled-up version of the formation of solar-type stars. 
Among the objects where rotation evidence has been inferred, 
one stands unique: this is G24.78+0.08 Al (hereafter G24 Al). 
G24 Al is one of the three massive rotating toroids detected in 
G24.78+0.08 (Beltran et al. [20041 |2005l >, a massive star forma- 
tion region located at a distance of 7.7 kpc. At the center of G24 
Al, an unresolved hypercompact (HC) Hn region has been de- 
tected at 1.3 cm by Codella et al. dl997l l. According to these 
authors, the spectral type of the star is at least 09.5, based on 
the free-free continuum flux, implying a stellar mass ~ 2OM . 
Very Large Array (VLA) observations in the NH3(2, 2) line have 
revealed that the gas in the toroid is undergoing infall towards 
the HC Hn region (Beltran et al. 2006), suggesting that accretion 
onto the star might still be going on, even through the HC Hn 
region as described in some models (see e.g. Keto & Wood 
2006). However, recent VLB A proper motion measurements of 
H2O masers associated with the HC Hn region (Moscadelli et 
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al. 120071 hereafter MGCBF) indicate that the latter might be ex- 
panding, thus questioning the possibility of accretion onto the 
star. 

To discriminate between the accretion and expansion scenar- 
ios, knowledge of the internal structure of the HC Hn region is 
crucial. According to MGCBF, the expansion should be driven 
by the stellar wind and the HC Hn region should thus have a 
shell-like appearance. With this in mind, we have imaged the 
free-free emission of the HC Hn region in G24 Al in order to 
resolve its structure. 



2. VLA Observations 

We carried out interferometric observations in the 1.3 cm and 
7 mm continuum with the VLA on April 7 and 8, 2006. The 
observations were carried out with the A array plus Pie Town 
configuration. The field phase center was set to the position 
a-(J2000) = 18 h 36 h 12 s :66, <5(J2000) = -07° 12' Iff.' 15. The detec- 
tion of strong H2O and CH3OH maser emission, respectively 
at 1.3 cm and 7 mm, within a few arcsec from the continuum 
source allowed us to self-calibrate on the maser spots and then 
apply the phase corrections to the continuum. The absolute posi- 
tion of the continuum at 1.3 cm was shifted 17 mas eastward 
and 35 mas northward to allow matching of our H2O maser 
spot positions with those measured with much greater accuracy 
by MGCBF with the VLBA. Final maps were created with the 
ROBUST parameter set equal to at 1.3 cm, and with natural 
weighting at 7 mm. The resolution and sensitivity (3<x) attained 
were 07095 x a/056 (PA=8°) and 0.25 mJy beam- 1 at 1.3 cm, 
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Table 1. Parameters of the 1.3 cm and 7 mm radio continuum of 
the HC Hn region in G24 Al 



Position" 










a(J2000) <5(J2000) 


A 




s r 


R b 


h m s '" ' " 


(cm) 


(K) 


(mJy) 


(AU) 


18 36 12.555 -07 12 10.83 


1.3 


4.8 x 10 J 


72 


580 


18 36 12.556 -07 12 10.80 


0.7 


2.4 x 10 3 


101 


550 



(a) obtained from a 2D gaussian fitting 

(b) geometrical mean of the major and minor semi-axes, obtained from 
gaussian deconvolution 



and 07061 x 07047 (PA=-4°) and 0.86 mJy beatrr 1 at 7 mm. 
Continuum images at 21 and 6 cm were obtained from VLA 
archive data (projects AB544 and AB515). The flux density at 
6 cm and the upper limit at 21 cm were calculated from gaus- 
sian fits to G24 Al and B (Codella et al. 1 19971 1. made with the 
task JMFIT of AIPS . Since the angular resolutions (57 1 x376 and 
575 x472 respectively at 6 and 21 cm) are insufficient to resolve 
Al from B (separation ~4"), in the fit we fixed the diameter of 
Al to the half power beam width (HPBW) of the corresponding 
map, and its position to that obtained from our VLA + Pie Town 
observations. All parameters of G24 B were left free. 

3. Results and discussion 

3. 1 . Structure of the HC Hn region 

Figure Q] shows the emission of the HC Hn region towards G24 
Al at 1.3 cm and 7 mm wavelengths. Table[T]reports the absolute 
position, the maximum brightness temperature measured in the 
synthesized beam, T s b, the flux density, S v , and the deconvolved 
radius, R, of the HC Hn region at both wavelengths. In Fig. [2] 
we show the radio continuum spectrum of G24 Al and the fit 
obtained for a classical Stromgren Hn region with radius of 07 1 3 
(1000 AU at 7.7 kpc) and a stellar Lyman continuum of 6.7 x 
10 47 s _1 , corresponding to an 09.5 type star. 

The HC Hn region appears to be almost circular at both 
wavelengths, despite a faint "tail" extending towards the south- 
west. This slight asymmetry could be due to lower density of the 
ambient gas, as suggested by the fact that the peak of the 1.3 mm 
continuum emission imaged by Beltran et al. (2004) lies on the 
opposite side of the HC Hn region (see cross in Fi.gQ]). Such a 
"tail" is unlikely to be related to the velocity gradient observed 
on a larger scale in the CH3CN lines (Beltran et al. 2004]!, as 
the directions of the two are quite different (PA^245° and 222° 
respectively for the "tail" and the gradient). Moreover, the ex- 
pansion of the HC Hn region cannot affect the velocity field over 
a region 10 times larger, like that where the gradient is observed. 

As can be seen in Figure Q] the emission is resolved at both 
wavelengths, and shows limb brightening, especially in the east- 
ern edge of the HC Hn region. The resolution at 7 mm is high 
enough to resolve the structure of the HC Hn region, which re- 
sembles a ring. To better study this ring morphology, a 7 mm 
map with ROBUST=0 and reconstructed with a circular beam 
of O/056 x O/056 was created. Using this map, we obtained pro- 
files of the emission at different angles by taking slices passing 
through the barycenter of the HC Hn region (cross in Fig. [3]), 
located at a(J2000) = 18 h 36 h 12 s :557, 5(J2000) = -07° 12' 10/80. 
Such emission profiles, not shown here, show two peaks, as one 
would expect in the case of an emission ring. In Fig. [3] we have 
outlined the ring by joining the positions of these peaks. This 
suggests that what we are observing in the continuum is an ion- 
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Fig. 1. Map of the VLA 1.3 cm (top panel) and 7 mm (bottom 
panel) continuum emission towards G24 Al. The contour lev- 
els are -3, 3, 6, 10, 15, 20, 25, 35, 60, 90, and 120 times a, 
where <x is 0.08 mJy beatrr 1 at 1.3 cm, and 0.29 mJy beatrr 1 at 
7 mm. The synthesized beams are drawn in the bottom left cor- 
ner. The black cross indicates the position of the 1 .4 mm contin- 
uum source (Beltran et al. 200"4il. 



ized gas shell. This is not in contradiction with the continuum 
spectrum of the region (Fig. [2]), because as discussed by Shull 
(1980), the free-free spectrum of a shell Hn region is basi- 
cally indistinguishable from that of a homogeneous, spherical 
Stromgren Hn region with the same angular radius and Lyman 
continuum. The only distinction between the shell and the spher- 
ical models is the presence of limb brightening at optically thin 
frequencies, as observed in G24 Al. 

Although the existence of an emission shell appears evident, 
the resolution of our observations is not sufficient to derive the 
intrinsic diameter of the ring in Fig. [3] For this purpose we have 
taken the brightness temperature profile along a cut (dashed line 
in Fig. [3} passing through the barycenter of the 7 mm image 
and the peak of the emission. This profile (see Fig. has been 
normalized both in temperature and radius, dividing the former 
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Fig. 2. Radio continuum spectrum of G24 Al. Values at 21 
and 6 cm are from VLA data archive, at 3.6 cm from Forster 
& Caswell ( 120001 1, and at 1.3 cm and 7 mm from this study. 
The triangle at 21 cm indicates an upper limit. The maximum 
brightness temperature measured in the synthesized beam is 
~ 5 x 10 3 K at 1.3 cm and ~ 2 x 10 3 K at 7 mm, indicating that 
the emission is rather optically thick. 

by the brightness temperature measured at the position of the 
barycenter (T B (0) - 1800 K) and the latter by the full width 
at half power (FWHP=0717) of the HC Hn region. For a shell 
Hn region, this profile depends on three parameters only: the 
ratio FWHP/HPBW, the optical depth r ff at the center of the Hn 
region, and the ratio RJR between the inner and outer radius of 
the shell. In our case, we know that FWHP/HPBW=3.085 and 
Tff = — ln(l — Tb(0)/T s ) =i 0.2 for an electronic temperature of 
10 4 K. The best fit is then obtained for Ri/R - 0.9 and is shown 
in Fig. a 

Such a result is in agreement with the prediction of 
Shull ( fl980T > that the shell should be very thin, i.e. (Ri/R - 1) <K 
1 . From the fit it is also possible to estimate the deconvolved ra- 
dius of the shell, which is 590 AU. Note that this is significantly 
smaller than the value of ~ 1000 AU needed to fit the observed 
radio continuum spectrum of the HC Hn region (see Fig. 
This is due to the fact that the HC Hn region is not perfectly 
spherically symmetric, but shows an emission "tail" towards the 
southwest (see above), which in spite of its relatively low bright- 
ness, conveys a non-negligible fraction of the total flux - ~28% 
(~20 mJy) at 1 .3 cm and -30% (~ 30 mJy) at 7 mm. This "extra" 
flux cannot be accounted for by the simple shell model. 

3.2. Expansion of the HC Hn region 

As shown before, the HC Hn region in G24 Al cannot be de- 
scribed in terms of a classical compact and spherical Hn region, 
but of an ionized shell. This result confirms the interpretation 
given by MGCBF to explain the distribution and proper motions 
of the 22.2 GHz water maser features towards G24 Al. As seen 
in Fig. [3] the water masers are well correlated with the eastern 
edge of the HC Hn region and have proper motions directed 
approximately perpendicular to the arc described by the maser 
spots. According to MGCBF, this configuration indicates that 
the water masers are tracing expansion, which in turn suggests 
that the HC Hn region is expanding. Because of the high velocity 
(~ 40 kms 1 ), such an expansion cannot be lead by the thermal 
pressure of the ionized gas, as in the classical case depicted, e.g., 
by De Pree et al. (Q995). Therefore, it must be some additional 




-1.3 -1.4 -1.5 -1.6 -1.7 

Aa (arcsec) 



Fig. 3. Map of the VLA 7 mm towards G24 Al {grey scale) 
created with ROBUST equal to and reconstructed with a cir- 
cular beam of 0'.'056 x 0'.'056. The contour levels range from 
1.3 to 9.1 in steps of 1.3 mJy/beam. The positions of 22.2 GHz 
water and 6.7 GHz methanol maser spots (from MGCBF) are 
denoted respectively by white and black circles. Black arrows 
show the measured absolute proper motions of the water maser 
spots (also from MGCBF). The white cross marks the position 
of the barycenter of the 7 mm continuum map. The white line 
connects the peaks of the profiles obtained from a large num- 
ber of slices through the barycenter (one every 10° in PA). The 
dashed line indicates the direction along which the normalized 
brightness temperature profile in Fig.|4]has been computed. 

mechanism driving the expansion. This is also suggested by the 
line width (40 km s _1 ) of the H66ff recombination line (Goddi, 
private communication): this is well in excess of the thermal 
width in Hn regions (typically 20 kms _I ) and thus requires ad- 
ditional broadening (see Sewilo et al. 2004), such as large scale 
motions, consistently with our expanding model. Note however, 
that the width of the recombination line appears to indicate a 
smaller expansion velocity (Robinson et al. 119821 Welty[l983) 
than that traced by the water masers. However, a caveat is in or- 
der. The problem posed by the detection of such broad recombi- 
nation lines is not yet solved (see Sewilo et al.|2004), especially 
because their strength is inconsistent with the Hn region being 
optically thick in the continuum. For this reason, to describe the 
recombination line emission from hypercompact Hn regions, 
and thus obtain a correct estimate of the expansion velocity, one 
needs a better model, where kinematics and radiative transfer are 
properly taken into account, which goes beyond the purposes of 
the present work. 

Following the model by Shull ( fT980b , MGCBF, conclude 
that the Hn region expansion may be driven by a powerful stel- 
lar wind. In such a model, there is an initial phase of free ex- 
pansion at the wind velocity, when the wind sweeps up its own 
mass in interstellar matter. After that, the system evolves into a 
four-zone structure, in which a thin, dense ionized shell contain- 
ing most of the swept up material is created. In this process two 
phases can be identified: a pressure-driven expansion, followed 
by a momentum-driven expansion after a critical time that de- 
pends on the wind mechanical luminosity and density of the sur- 



4 Beltran et al.: The hyperyoun; 



FWHP/HPBW=3.085 T ff =0.2 R/R s =0.9 
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Fig. 4. Normalized brightness temperature profile (histogram) 
along the dashed line in Fig. [4] The offset is normalized with 
respect to the measured FWHP. The curve represents the best fit 
obtained with a shell model of the Hn region - see text. 

rounding environment. For a given wind mechanical luminosity 
and density, the radius (Run) and velocity ( Vhii) of the expanding 
shell can be expressed as a function of time (Shull 1980 ). For our 
estimates we have adopted a Lyman continuum of 6.7 x 10 47 s _1 
(see Sect. 13.11 ). an inner-to-outer radius ratio of 0.9 for the shell 
(see Sect. 13. U . and the fiducial values of 10 7 cirT 3 , 10 36 erg s" 1 , 
and 2000 km s respectively for the environment density, me- 
chanical luminosity of the wind, and wind speed (see MGCBF 
for a detailed discussion on the input model parameters). 

The values of the shell radius and velocity as a function of 
time are shown in Fig. [5] as well as the 1.3 cm and 7 mm con- 
tinuum fluxes of the HC Hn region. For the sake of comparison, 
we also show the same quantities for a classical Stromgren Hn 
region (dashed curves). The corresponding observed values with 
the associated uncertainties are also reported and demonstrate 
that only the wind-driven can fit all the data. Note that the fluxes 
used here (56 mJy at 1 .3 cm and 74 mJy at 7 mm) do not contain 
the contribution from the emission "tail" mentioned in Sect. 13. II 

Given the very short time scale involved (21-66 yr), the size 
variation expected in only 5 yr (~10 mas) should be easy to re- 
veal using the H2O maser spots as probes, thus providing us with 
a robust test for our model. 

4. Conclusions 

Our radio continuum images with <90 mas resolution have re- 
vealed a ring shaped structure in the HC Hn region at the center 
of the rotating toroid imaged by Beltran et al. ( 120041) . supporting 
the model of a wind-driven, ionized shell proposed by MGCBF. 
We thus conclude that very likely the HC Hn region is expanding 
on very short time scales (21-66 yr). This finding appears to con- 
tradict the fact that the gas on a larger scale is infalling towards 
the HC Hn region (Beltran et al. 2006). A possible explanation 
is that the infalling gas is not accreting any more onto the star, 
but is stopped (and possibly diverted into the outflow observed 
by Furuya et al. 120021 ) at the surface of the HC Hn region, right 
at the shock front traced by the H2O maser spots. It is worth not- 
ing that the CH3OH masers observed by MGCBF appear to lie 
further from the HC Hn region than the H2O masers, suggesting 
that they might be located in the pre-shock material. If this is the 
case, methanol masers should be participating to the infall and 
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Fig. 5. Flux (top panel), velocity (middle panel), and radius 
(bottom panel) of a HC Hn region as a function of time for 
two different models: a classical Hn region (dashed lines), and a 
wind-driven shell (solid lines). The thin lines in top panel indi- 
cate the flux at 1.3 cm, and the thick lines the flux at 7 mm. The 
measured fluxes (56 mJy at 1.3 cm and 74 mJy at 7 mm) do not 
take into account the contribution from the southwestern emis- 
sion "tail". The vertical error bars correspond to uncertainties of 
30% for the radius and expansion velocity, and 15% for the flux. 
The error bar on t corresponds to the uncertainty estimated by 
MGCBF. 

their proper motions should be directed towards the HC Hn re- 
gion, unlike the water masers which are clearly expanding from 
it. Future VLBI multi-epoch measurements are bound to confirm 
or rule out this scenario. 
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